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SECTION I 


CHI SQBARE D PROGRAM 


I. 


Introduction 

1.1 Purpose 

To evaluate CHI SQUARE. 

1.2 Method 



m 

2 

k=l 


(Ok - Ek>^ 

Ek 


Where Oj^ is the observed frequency. 

Ejj Is the expected frequency, 
m is the number of categories. 
II. Usage as a Routine 

2.1 Minimum System Components Required 
Fortran Compiler D System, 


2.2 Input 

The Input data Is read from card reader. The first card of 
each problem must be a control card, followed by data cards, 
2,2,1 First Control Card Format 


Coltunns 

Contents 

Description or Remark 

1-5 

Blank 


6-10 

Problem ID, 

An Integer < 99999 

11-20 

No. of categories 


21-80 

Not used 



1-1 


2 . 2.2 


2,2.3 


2.2.4 


Data Card Format 


Columns 

Contents 

Description or Remark 

1-5 

Blank 


6-10 

Input Data ID 

This number should be 
exactly the same as 
the number vihich appears 
In column 6-10 in the 
control card. 

11-20 

Category Number 

An integer. 

21-30 

Observed frequency 
for kth category. 

Decimal 

31-40 

Expected frequency 
for kth category. 

Decimal 

41-80 

Not used 


Last Control Card Format (Thi 

the 

s card is needed only at 
end. of a set of problems) 

Columns 

Contents 

Description or Remark 

1-5 

Blank 


6-10 

99999 


11-80 

Not used 



Card Sequence Layout 


2nd ProhLen 
data 


First Problem 
data 
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2,3 Output 




The output Is printed on-line, 
III, Example 

3,1 Problem No, 1 


Category 

Observed Frequency 

Expected Frequency 

1 

90 

100 

2 

53 

50 

3 

4o 

30 

4 

45 

50 

5 

12 

10 


3.2 Input Card Format 



OOOOOOOOOOOOOOQOOOOQOOOQOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOODOOOOI'OOOOOOODODOODOOOOOO 

<73 4 5 6 7 8 6 15 11 12 13 14 15 16 171619 70 712273 24 75 26 27 7673 30 31 37 33 34 3536 37 38 36 40 61 47 4344 45 46 47 48 43 50 51 57 53 51 55 56 57 56 59 60 61 62 63 64 65 66 67 68 66 70 71 72 73 74 75 76 77 76 79 10 

1 11 1 11 1 n|t 1 1 11 1 1 1 1 1 11 1 1 1 n 1 i 1 1 1 < 1 1 1 1 1 1 1 ! 1 11 M 1 1 11 1 1 n 1 11 11 1 1 1 11 1 11 1 1 1 1 11 1 1 1 1 in 
2222222222222222222222222222222222222222222222222222222222222222222222222222222? 
33333333333333333333333333333333333333333333233333333333333333333333333333333333 
444444444444444444444444444444444444444444444444444444444444444444«4444444444444 
3 S 5 SS 5 S 55 SS 5 SSSSSSS| 5 S 9 SSSSSS 55 S 55 S 5 S 95 S 5 S 35 S 555 S 5 S 555535555 S 553555555 S 5 S 3 S 55553 
66 E 66 S 66 G 66666666666666 S 666666 S 6666666666 G 66666 S 66 E 6668666 C 6666 E 66 E 6 E 666 B 6 S 66666 

7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 n 7 7 in 7 7 7 ? 7 7 7 7 7 7 7 7 7 7 7 7 7 7 ? 7 7 7 7 7 7 7 7 7 7 7 

8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 S 8 8 8 G 8 8 8 8 8 8 8 S 8 S 8 8 8 8 E S 8 8 S 8 8 S 8 8 8 8 8 8 8 8 8 8 8 8 8 8 e 8 8 8 8 8 8 8 8 8 8 8 8 8 

9 9 9 S 9 9 9 9 3 9 9 9 9 9 9 9 9 9 9 9 9 9 9 S 9 9 9 9 9 9 9 S 9 9 9 9 S 9 9 9 9 ? 1 9 9 9 9 3 3 5 9 9 S 9 S 9 9 S 9 P S 9 9 9 9 ? 9 9 9 3 9 9 9 9 9 9 S 9 9 9 

• ? 3 < 3 6 7 8 9 10 II 1? 13 14 15 16 17 18 19 20 21 2; 23 24 25 26 27 28 29 3C 31 32 33 3* W3«3T 38 3fl40*'<?43«85<6 4’ «3«5D 57 M55M5?53 5t‘l)1) 61 « 63 5*65 BES7 68 89 70 71 U 73 74 71>rK T7 7g?ci 8p 


3 . 3 Output 

CHI SQUARE PROGRAM 
PROBLEM NO. 1 
NO. OP CATEGORIES » 5 
CHI SQUARE » 5.^1333 
ALL DONE. 


IV. 


Disclaimer 

This program has been checked hut no guarantee Is expressed or 
Implied to the user concerning Its functioning. Honeyirell will 
cooperate In making functional those programs written by them. 
However, they assme no responsibility for programs distributed 
but not written by Honeywell. 
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SECTION II 


LEAST SQUARES CURVE FITTING D PROGRAM 


I . Purpo se 

Fits a polynomial of degree n to a set of m observations by the 
method of least squares. It Is written In Fortran language, 

II. Method 

2.1 Symbol 

n = degree of polynomial 
m = number of observations 

Xj = J-th observation of independent variable 
yj = J-th observation of dependent variable 

A 

yj = predicted value of dependent variable 

^3 = deviation between observed and predicted value of 
J-th observation of dependent variables, i.e, 

Xj = 1-th power of Xj 

2.2 Description 

For a given n observations of (xj, y.) to fit a polynomial 
of the following form: 

y = bg + b^x + bgx^ + ... + b^x’^ (1) 

This program Is used to determine the coefficients bi*s 
by the method of least squares. 

The principle of Least Squares seeks to minimize the squares 
of the residuals „ m n 

ZS = S (y. - E BiX^)^ 

J=1 1=0 ^ 


2-1 


where is the estimated value of in (1). 
The normal equations are: 

hgExj + + . . . + 

b(,rxj + + ... + 


= Sxjyj 


• • • • ( 2 ) 

• • • • 

• • • • 

This is a set of (n+1) simultaneous linear equations in 
hj^ which is solved by Gaussian Elimination in this program, 

2,3 Matrix Operations 


The augmented matrix which is the coefficient matrix of (2) 
with the right hand side as an additional column is denoted 
as A = (aij). All (n+l)(n+2) elements will be stored in 
memory. 


The diagonal element is used as a pivot. Denote the pivot 
as ajf The new matrix will be generated by the following 
algorithm: 



for l=i=k 

(3) 


for l=k, ^?^k 

(4) 

^IJ “ ^lk^kj/®-kk 

for lA, J?^k 

(5) 

-aik/ajjk 

for i/k, J=k 

(6) 


where a^j is the element of the new matrix, 

2,4 Estimates of bj^ 

After n+1 iterations of the matrix operation in 2,3, the 
last column gives the estimates 
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2,5 Calculation of predicted values and deviations 


III. 


IV. 


V. 


VI. 


The predicted value for the dependent variable Is 


A n 

y. = 2 E 
1=0 



The deviation between the actual and predicted values of 
the dependent variable is 



Minimum System Components Required 


3.1 Fortran Compiler D 


3.2 One work tape 


Input 

The input data Is read from the card reader. The format 
will be given In detail In section 8. 

Output 

All output is printed on-line, A sample output printout 
will be shown in 9.3, 

Estimate of Maximum Degree 

(N+l)(N+2)(K+2) < M - 13,600 

where N = degree of polynomial 

K = length of mantissa of floating point numbers 

M = memory configuration 

As an example: Assume K = 10 K N 

T&k “iT 

20K 18 

24k 27 

28K 32 

32K 37 

This version Is written for 16K with N=12, precision 10. In 
order to adapt this to another memory size, the dimension of V 
must-be changed in the “common” statement of each chain. 
Example: For 32K, precision 10, the dimension of V should be 

V(N+1, N+2), l.e., V(39, 40). 
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VII, Preparation of Input Data 


7.1 Size of Independent variable 


The smallest value of the independent variable Xj must 
satisfy the following condition 

> 10 "^^ 

For example, to fit a polynomial of degree 16, the smallest 
value Xj must be greater than 10”>’, roughly. 


The largest value of the independent variable xj must 
satisfy the condition 


S X 2n 

J=1 ^ 


^ 10 


99 


For example, for a given 100 observations of (xj, y^ ) to 
fit a polynomial of degree 16, the largest value xj should 
be less than 10 3, roughly. 


7,2 Number of problem 


A maximum of 99998 problems can be run at one time. 


7.3 Nmber of observations 


The number of observations that can be handled is one less 
than 10^ and subject to the condition that 

S x,2n < io99 

J=1 ^ 

which is mentioned In 8,2, 

7.4 Card Format 

8,4,1 First control card of each problem 


Columns 

Contents 

Description or Remark 

1-5 

Blanks 



6-10 

Problem ID, 

An 

Integer <99999 

11-15 

Degree of Polyno- 
mial 

An 

integer, right 

16-20 



Justified 

Total Number of 
Observations 

An 

Integer 
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Colximns 

Contents 

Description or Remark 

21-24 

Blanks 


25 

Blank or 1 

If "I", program will 
suppress matrix 
elements printout. 

26-80 

Not used 


Data card 



Columns 

Contents 

Description or Remark 

1-5 

Blank 


6-10 

Input Data ID 

This no. should be 
exactly the same as 
the no, which appears 
in Column 6-10 in the 
Control Card. 

11-20 

Observation No. 

An Integer, 1 stands 
for first observation, 
2 for second observa- 
tion, etc. Order is 
irrelevant. 

21-30 


Decimal, value of 
Independent variable. 
The decimal point is 
assumed at extreme 
right, if there is no 
decimal point present. 

31-40 


Decimal, value of 
dependent variable. 
Decimal point is the 
same as Xj field. 

41-80 

Not used 

Last control card 


Coliimns 

Contents 

Description or Remark 

1-5 

Blank 


6-10 

99999 

An Integer. 

11-80 

Not used 
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7.^.4 Card Sequence Layout 


2nd Problem 
data 


1st Problem 
data 


last control 
card 




control card 
of problem No. 2 


i> * ' 
'‘‘1 
Ir 


Second data card 


First data card 
of problem No.l 




control card 
of problem No*.l 


VIII. Sample Problem 

8.1 Find by the method of least squares a formula of the type 

y = bQ + b^^x + b2X^ + b^x^ 
which will fit the following data (homogeneous) 


X 

-4 

-2 

-1 0 

1 

3 

4 6 

y 

- 35.1 

15.1 

15.9 8.9 

.1 

.1 

21.1 135.0 


8.2 


Input Bata 



OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOaOOOOQOOOIlOOIlDOOOOOOOOOOOOOOOOOOOOOOOOOOO 

1 2 3 4 5 8 7 I 9 10 II 12 13 14 15 IS 17 I1 18 20 21 22 23 24 25 28 27 2829 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 48 47 48 4S 50 51 52 S3 54 5S 56 57 51 59 80 81 82 83 64 65 68 67 88 89 70 71 72 73 74 75 78 77 78 79 80 

n 1 11 1 n 1 11 1 1 1 1 11 1 11 1 1 1 11 11 1 n 11 1 n n 1 n 1 1 1t 1 1 1 11 11 n 11 1 1 1 111 11 1 11 1 1 1 1 n 1 1 1 1 11 1 1 


22222222222222222222222222222222222222222222222222222222222222222222222222222222 
33333333333333133333333333333333333333333333333333333333333333333333333333333333 
44444444444444444444444444444444444444444444444444444444444444444444444444444444 
55S555553IS555S595555555S5S5S5S555S55555555555555S5555S5555555SSS5S5S55555555555 
S666E666E66666E6EE666S6666666686SE6686666G66EE6E666e66EEEEE666BES6EB6EE6EEEE6EEE 
7 7 7 7 7 7 J 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 J 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 J 7 7 7 7 
1 8 8 8 8 8 S 8 8 S S E 8 S 8 S 8 S 8|S 8 8 8 8 3 8 3 E 8 8 8 8 8 a 8 8 8 8 S 8 8 8 8 8 S 8 8 8 3 S 8 8 B 8 8 8 8 8 8 8 S 8 8 8 8 S 8 S S B 8 S 8 8 8 8 8 8 8 
99939999999999999S9999999999999993999399993999993S999999999999999999999999999999 

12 3 4 5 6 7 8 8 10 11 12 13 14 15 18 17 14 15 20 21 22 23 24 25 26 27 26 39 30 3U2 33 34 35 36 37 38 39 40 41 42 43 44 45 48 47 4845 50 51 52 53 54 55 58 57 58 58 SO 81 82 63 64 85 68 67 68 88 TO 71 72 73 74 75 78 U 76 78 80 


2-7 


LEAST squares CURVE FITTING PROGRAM (LESCUF) 


PROBLEM NOf 5 

NO, OF observations = 6 

DEGREE OF POLYNOMIAL * 3 


ELEMENTS OF A-MATRiX 


A ( 

1« 

1) 

E 

0.800000000E 

01 

A( 

If 

2) 

s 

0.700000000E 

01 

A( 

If 

3) = 

0.830000000E 

02 

A{ 

If 

4) 

s 

0.235000000E 

03 

A ( 

If 

5) 

X 

Of 161100000E 

03 

A ( 

2f 

2) = 

0.830000000E 

02 

A« 

2t 

3) 

s 

0.235000000E 

03 

A( 

2f 

4) 

X 

0.190700000E 

OA 

A ( 

2f 

5) 5 

0.989100000E 

03 

A ( 

3 f 

3) 

r 

0,190700000E 

04 

A ( 

3f 

4) 

X 

0.798700000E 

04 

A ( 

3f 

5) s 

0,471330000E 

04 

A( 

4f 

4} 

s 

0.556430000E 

05 

A( 

4f 

5) 

= 

0»326229000E 

05 

A( 

f ) 





CONSTANT TERM = 0,901104386E 01 

POWER COEFFICIENT 


to 

I 

00 


1 

-0f8966l43l9E 

01 



2 

-Of 100009408E 

01 



3 

0.999074298E 

00 





PREDICTED VS 

ACTUAL RESULTS 


RUN NC, 

VALUE OF X 

actual 

predicted 

DEVIATION 

1 

•4,000000 

-35,100000 

-35,066644 

-,033356 

2 

-2,000000 

15,100000 

14,950360 

,149640 

3 

-1,000000 

15,900000 

15,978019 

-,078019 

4 

,000000 

8,900000 

9,011044 

-,111044 

5 

1,000000 

,100000 

.043861 

,056119 

6 

3,000000 

,100000 

,086774 

,013226 

7 

4,000000 

21.100000 

21,085721 

.014279 

8 

6,000000 

135.000000 

135,010846 

-,010846 

END OF CALCULATIONf PROBLEM NO 

. 5 




ALL DONE, thank YOU, 


8.3 Output Data 
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SECTION III 


MEAN, VARIANCE AND CORRELATION D PROGRajl 


I . Purpo se 

This program is used for computing mean, variance, covariance, 
standard deviation and correlation coefficient of the variables 
which are stored in two groups on the tape, the X-group and 
Y-group. 

II. Method 

2.1 Symbol 

m = number of observations 
p = nmber of variables in X-group 
q = number of variables in Y-group 
= i-th observation of variable xj 
yij ~ i-th observation of variable yj 
= 1-th observation of variable zj 
Xj = Zj (for 1 S J ^ p) J-th variable in X-group 
yk = (for 1 < k s q) k-th variable in Y-group 

Zj = mean of variable Zj 
<Tzj = variance of Zj 
CTzj = standard deviation of Zj 
cov(xj, yjj) = covariance of xj and y^^ 

^k = coefficient of correlation between Xj and y^^ 

2.2 Formula 

For each variable, we compute 

Zj = Z Zj^j/m 
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(Tzj = (E - razj^)/(m-l) 



For each pair of variables, l.e, one of which Is In the 
X-group and one In Y-group, compute 

cov(xj, = (E - mXj yi^)/(m-l) 

and . 

ojk = cov(xj, yif)/(^j 


III, Minimum System Components Required 
3.1 Fortran Compiler D System. 


IV, Input 

The input data is read from the card reader. The input data 
format villl be given in detail in section 8.6, 

V, Output 

All output is printed on-line, A sample output will be shown 
in section 9» 

VI. Subroutine 

The square root routine must be supplied from the library of 
subroutines. 
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High Speed Memory Requirement 



Memory Space 

Notation Used 
in Program 

Remarks 

A 

p • q 

VECTOR(I,J) 

matrix A space has multiple 
uses as: 

1) sum of cross products 

2) covariance of variables 

3) correlation of coeffi- 
cients 


p 

Z(I) for l<I<p 

Z(I) space has multiple 
uses as: 

1) inputs x^, y^ 

2) residiaal sum of squares 

^i 

q 

Z(l) for p<l<I<p+q 


fji 

p 

SUM Z(I) for l<I<p 



q 

SUM Z(I) for p+l^I<p+q 

ZY ^ 
J ^ 

p 

SUM ZZ(I) for 1<I^ 


q 

SUM ZZ(I) for p+l<I<p+q 

<r± 

p 

SIca(A(l) for ksl<ip 



q 

for P+l<I<p+q 


"i 

p 

Z BAR(I) for 


«• 

q 

for p+l<I^+q 



VIII. Preparation of Input Data 


8.1 Size of data 

If the largest sum of squares or the sum of the cross 
products is greater than 10^, the printout xflll lose the 
most significant digits. However, the statistical results 
will be correct. The size of data should depend on the 
number of observations, 

8.2 Number of variables 

[p*q + 5(P+q)D (d+2) < M - 14,000 
where M is the memory size 

p is the number of variable in X-group 
q is the number of variable in Y-group 
d is the number of digits accuracy. 

As an example: p = q d = 10 

Memory Size Estimated Maximum Value For p 


16K 


10 


24k 25 

32K 3^ 


* Remark: This version is written for 16K machine with 

precision 10, In order to adapt to other size 
machines, the arrays in the "common'* statement 
of each chain must be changed according to the 
above estimated formula, 

8,3 Number of problems 

A msoiimum of 99998 problems can be run at one time. 
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8.4 Number of data points 

The number of data points (observations) that can be handled 
is less than 10^ and is limited by the restriction that the 
largest sum of squares must be less than 10^. This program 
handles an equal number of data points for each variable. 

8.5 Decimal point 

The decimal point of input data need not be entered if the 
number is an integer. The decimal point will be assumed 
to be at the extreme right. For example, columns 21 to 30 
are used for the first variable, say xi, within a data card. 

If a **15” is punched in columns 29 and 30, then x;j^ is 15 • 

If **15” is punched in columns 28 and 29, then x^ will be 
150 . A number with a decimal point can be punched anywhere 
within the specified colmns for that variable. For example 
a ** 25 •4** in columns 19 to 23 is the same as if punched in 

26 to 30 . 

8.6 Tape format 

8.6.1 First control card of each problem 

Columns Contents Description or Remark 

1-5 Blank 

6-10 Problem Identification An integer <99999 

11-15 Nmber of observations An Integer 

16-20 Number of variables in An integer 

X- said Y-groups. 

21-25 Number of variables in An integer * If 

Y-group. this no. is zero, 

the program will 

not calculate corre- 
lation. 

26-80 Not used. 
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8.6.2 Data card 


Columns Contents 
1-5 Blank 

6-10 Input Data ID 

11-15 Logical Sequence No. 


16-20 Observation No. 


21-30 

31-40 zg 

41-50 Z 3 

51-60 

61-70 Z3 

71-80 


Description or Remark 


This no. should be exactly 
the same as the no. which 
appears in Col. 6-10 in 
the Control Card. 

An integer. 1 stands for 
first 5 variables data 
card, 2 for second 5 
variables data card, etc., 
if the nmber of variables 
is less than 5» a 1 is 
sufficient for all data 
cards. 


An integer. 1 stands for 
first observation, 2 for 
second observation, etc. 



Decimal, value of first 
data within a card 

Decimal, value of second 
data within a card 

Decimal, value of third 
data within a card 

Decimal, value of fourth 
data within a card 


Decimal, value of fifth 
data within a card 

Not used 
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8.6.3 Last Control Card 



Columns Contents Description or Remark 

1-5 Blanks 

6-10 The Integer 99999 

11-80 Not used 


8.6.4 Card Sequence Layout 
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IX. Sample Problem 
9 . 1 Problem 2 


Given the following data, compute the mean and variance 
of each variable and covariance and correlation coefficient 
between the variables in X-group and Y-group. 


m No, 

^1 

X, 

c 

1 

7. 

26. 

2 

1. 

29. 

3 

11. 

56. 

4 

11. 

31. 

5 

7. 

52. 

6 

11. 

55. 

7 

3. 

71. 

8 

1. 

31. 

9 

2. 

54. 

10 

21. 

47. 

11 

1. 

4o, 

12 

11. 

66. 

13 

10. 

68. 


X3 

^1 

Y2 

6. 

60 . 

78.5 

15. 

52. 

74.5 

8. 

20. 

104.3 

8. 

47. 

87.6 

6. 

33. 

95.9 

9. 

22. 

109.2 

17. 

6. 

102.7 

22. 

44. 

72.5 

18. 

22. 

93.1 

4. 

26. 

115.9 

23. 

34. 

83.8 

9. 

12. 

113.3 

8. 

12. 

109.4 
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9# 2 Input Data Cards For Card Image 



0 fi 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 D 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 9 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 

1 2 3 4 i t I t 9 10 111213 >4 B« nm»20» 222324ff2«2T»29 303U2»M3S3t 37»3S40 41 42 4344 4S4l4;4t4tS0SI S2S3S4SSS(S7SISt(0 t162(3MU»67S8ll7071 n» 74 7S7t n 717110 

11 11 1 1 1 1 1 1 n 1 11 ti i|ii 1 11 1 1 1 11 1 11 1 11 1 1 1 11 1 1 1 1 1 1 1 1 1 1 1 11 11 1 1 1 1 11 11 n 1 11 1 1 1t 1 1 11 1 1 11 
222222222|22222222222222222222222222222|2222222222222222222222222227222222222222 
33333333333333333331333333333333333333333333333333333333333333333333333333333333 
44444444444444444444444444444444444444444444444444444444444444444444444444444444 
5 S 5 5 5 5 5 5 5 5 5 S 5 5 5 5 5 5 5 5 S 5 5 S 5 5 5 5 5 1 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 S 5 5 5 5 5 5 S 5 5 5 5 5 S 5 5 5 5 5 S 5 5 5 5 5 5 5 5 5 5 5 5 
686S66666S6S66 6666666 B6 BBS 66SE6 66 B66S6666E68B&666 66 66 B6S66 66 666666 66 6666 66 66666 E 
77777777777777777777777777777777777777777777777777777777777777777777777777777777 
8888888886868888BB888686B86i8BSB888888S8888888888B888888B888i8!S!88888888SS88868 
99999939993999999999998989999$9999999399999999939939999999993il999S99S999999S9399 

I ] 3 4 i { 7 I llOI1t3UHB»l7l*ll»n2223»K»n»:9 30 31 U U 34 3S3« 37 3t 39 40 41 43 43 44 45 46 47 4141 N $1 S2 S3 54 55 56 57 51 31 60 61 67 63 84 6S 66 67 86 d TO 71 73 73 74 75 76 77 71 71 60 









9,3 Output Data 

• ••STATISTICS program - MEAN, VARIANCE AND CORREL AT I ON«»» 


PROBLEM NO, 2 

NO, OF observations = 13 

NO, OF variables in X-GROUP s 3 

NO, OF variables IN Y-GRQUP e 2 


MEAN 

7,A61536 

46,153646 

11,769231 

30,000000 

95,423077 


VARIABLE SUM SUM OF SQUARES 


X( 

1) 

97,000000 

1139.000000 

X( 

2) 

626,000000 

33050,000000 

X ( 

3) 

153,000000 

2293,000000 

Y( 

1) 

390.000000 

15062,000000 

Y( 

2) 

1240,499999 

121088,089997 


SUM OF CROSS products 


X ( 

1) VS Y( 

1) • 

2620,000000 

X ( 

3) VS Y( 

1) ■ 

4628.000000 

X ( 

2) VS Y{ 

2) » 

62027,799999 


VARIABLE 


VARIANCE 


X ( 

1) 


34,602564 


X ( 

2) 


242.141026 


X ( 

3) 


41.025641 


Y( 

1) 


280,166667 


Y ( 

2) 


226.313590 








COVARIANCES 

X( 

1) 

VS 

Y( 

1) 

a 

•24,166667 

X( 

3) 

vs 

Y ( 

1) 

B 

3,166667 

X ( 

2) 

vs 

Y ( 

2) 

s 

191,079467 







CORRELATION 

X( 

1) 

vs 

Y ( 

1) 

M 

-.245445 

X( 

3) 

vs 

Y( 

1) 

m 

,029537 

X( 

2) 

vs 

Y ( 

2) 

s 

,816253 


X I 

2) 

VS 

Y ( 

1) * 

15739,000000 

X ( 

1) 

VS 

Y ( 

2) ■ 

10031,999999 

X ( 

3) 

VS 

Y ( 

2) « 

13981,499999 


standard deviation 

5,682394 

15,560881 

6,405126 

16,738180 

15,043723 


OF variables 


X ( 

2) 

VS 

Y( 

1) 

S 

•253.416667 

X ( 

1) 

VS 

Y( 

2) 

B 

64.663462 

X ( 

3) 

vs 

Y ( 

2) 

m 

-51.519231 

COEFFICIENTS 





X ( 

2) 

vs 

Y ( 

1) 

m 

-.972955 

X ( 

1) 

vs 

Y( 

2) 

8 

.730717 

X ( 

3) 

vs 

Y ( 

2) 

8 

-.534671 


END OF calculation, PROBLEM NO, 2 


all done, thank you. 
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9*4 Special Comments on Error Conditions 

If any or all of the following comments appear on output 

9.4.1 DATA ID DOES NOT AGREE WITH CONTROL CARD. INID = 
xxxxx xxxxx NOC = xxxxx INNO = xxxxx 

9.4.2 DATA CARD OUT OF LOGICAL SEQUENCE. CHECK. INID = 
xxxxx xxxxx NOC = xxxxx INNO = xxxxx 

9.4.3 data card out of logical SEQUENCE. INID = xxxxx xxxxx 
NOC = xxxxx INNO = xxxxx 

then these errors are caused by the mispunching of the 
identiflc&tipn number of test data of a problem or the 
misplacing of data cards. The program will continue but 
the user should recheck the Input data deck. The numbers, 
INID, NOC, and INNO will help the program to pinpoint the 
mlspunched card or misplaced card. After correcting the Input 
data, the program should be re-run. 

Timing 

There are (m+l)(p+q) + pq additions, (m+3)(p+q) + 3P9 multiplica- 
tions, pq + 3 divisions and p + q square root subroutine calls. 


Disclaimer 

This program has been checked but no guarantee Is expressed or 
Implied to the user concerning Its functioning. Honeywell will 
cooperate In making functional those programs written by them. 
However, they assxime no responsibility for programs distributed 
but not written by Honeywell. 
















CT- 















SECTION IV 


STEP-VrtSE MULTIPLE REGRESSION ANALYSIS D PROGRAM 


I. Purpose 

This program is used to find the best fit of an equation of the 
following form: 


y = Bq + + BgXg + ... + (1) 

where y is the dependent variable and X 2 , ..., are independent 
variables. 

II. Method* 

2.1 Symbols 

n = number of Independent variables and dependent variables 
m = nimiber of sets of observations 
FI = F value for entering the variable 
F2 = F value for removing the variable 
TOL = Value of tolerance, an arbitrary constant 
Xijj = k-th observation of i-th variable 

^hk = yk = k-th observation of dependent variable 
Xj_ = mean of 1-th variable 
x^ = mean of dependent variable 
sij = residual sum of squares and cross products of 1-th 
and J-th variable. 

<Ti = \/s^ 

2TlJ = simple correlation coefficient of 1-th and J-th variable 
BjL = true value of coefficient of 1-th variable 
bQ = constant term of regression equation 
= estimated coefficient of 1-th variable 
Sy = standard error of dependent variable 
st)i = standard error of coefficient of 1-th variable 
^k = predicted value of dependent variable for k-th obser- 
vation. 

= deviation between observed and predicted value of k-th 
observation of dependent variable, l.e. y^^ - 
Njjjin, = subscripts of selected independent variables 


* Ralston, A. and Wllf, H. S., Mathematical Method for Digital 

Computers, John Wiley, New York, 

i 960 , pp. 191-203 
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N » subscripts of an Independent variable Just removed from 
regression equation 

^min ■ variance increase by deleting a variable 

“ variance reduction by adding a variable N 

max max 

^ » the degree of freedom 
2.2 Description: 

Let y be an estimate of y; the estimating equation may then be 
expressed as: 


r - »o ♦ *>1 \ * »2 ^ * Vl Vl 

where the b's are to be determined by means of available data. 

For a given m sets of values of n variables^ it is always possible 
to obtain b's^ the estimates of the coefficients B's. Geometrically, the 
problem is one of finding the equation of a hyperplane which best fits 
a set of m points in n dimensions. The problem now is to find the set of 


b's in (2) that will minimize the sum ^ * 

k-1 k 

In the step wise procedure, we obtain the regression equation by 

adding one variable 4t a time. The Intermediate iregression equations are: 

A 


■ b ' + b ' 

0 1 

X 

kl 

(3) 

- b" + b" 
0 1 


(U) 


* 

4 


The variable added is that one which makes the greatest inprovement in 
"goodness of fit". The coefficients are the best values when the regression 
equation is fitted by the specific variables which have already been entered 
in the equation. A variable, s^y x^^* ^ significant variable in an 
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early stage and thus enter the equation. But after several other variables 
have entered the regression equation, the “®y ^ insignificant. This 
variable will be removed from the regression equation before entering 
an additional variable. The final regression equation includes only those 
variables which are significant. 

The b's with the superscripts are intermediate estimated values of the 
coefficients of variables which are in the regression equation at that stage. 
There is no relation between the b's with the lower and the higher superscripts. 
The b of equation (2) will be equal to the b with the superscript of the 
coefficient of Xg in the final regression equation. And = 0 for the 
corresponding 1-th variable is not in the final regression equation. For 
example, we e^qpect to have the estimate of y: 

y * h * h * *>2 *2 * h *3 * h ^ * N ^5 

for n ■ 6 

If we have the final regression equation of the following foz*m: 

y»b "'+b. "'x- + b^"' Xi +b,"' x^ 

^ o I22u35 

then bg - bo"' , bg - b^"' , b|^ - bg*" , b^ - b^"' 

and b^ ■ b^ ■ 0. 

2.3 Basic Calculations and Formulas 

2.3.1 Calculation of sums of squares and cross products 

M 

Sums of variables for i ■ 1, 2, ... n 

B 

Sums of squares and cross products Xj^ for i«l, 2, ... n 

k»l j»l, 2, ... n 
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2*3.2 Calculation of means and residual sums of squares and cross products 

m 

Means % ■ ^ x /m for i«l, 2, ... n (5) 

k»l ik 

Residual svira of squares and cross products 
m mm 

8, . - ^ X 5~ X / ■ 


k«l ik Jk k»l ik k-1 Jk 


n - m 

■ 5~~ XX - X X for i-1, 2, ... n (6) 

k-1 ik jk i k-1 Jk j-1, 2, ... n 


2.3.3 Calculation of partisQ. correlation coefficients 

-\Js^ for i - 1, 2, ... n (7) 

y' ■ S ! (T^ ^ 

^ ij ij i ^ j for i - 1, 2, ... n, j«l,2,... n (8) 

y ■ Y y • 1.0 for i-1, 2, ... n (9) 

^ ji ij "ii 


2.3.U 


Normal Equations 
By minimizing the sum 



we have the normal equations 


ij J ^ for i-1, 2, ..., n-1 (10) 

They comprise a set of n-1 simultaneous linear equations in b which 
can be solved by any method. In this program we adopt the Gaussian 
elimination method. In the stepwise procedure, we obtain the final 
regression by adding one variable at a time. How to enter or remove 
a variable from the regression will be discussed in 2.3*5»s* 

2*3*5 Calculation of new matrix 

The original matrix consists of coefficients of correlation, 
i.e. k • }. It is symmetric. All n^ elements will be stored 

in memory. 
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Calculation of 

can be computed from matrix elements as follows: 

V . ■ a a /a^^ for each a, . ^ TOL 

i in ni ii ii ^ 

Where a^^ Indicates an element of Matrix A* 

A tolerance is placed on a^^ since there is a possibility of 

degeneracy in the matrix when an independent variable is approximately 

a linear combination of other independent variables. A variable is 

not allowed to enter the regression unless the a associated with 

ii 

it is greatet* than a preassigned value. The value of TOL is generally 
assigned a value which lies between .001 and .0001. In this program 
the value is assigned by the user in the control card. See 8.6.1. 

Another reason for establishing a limit for a^^ is to 
prevent arithmetic overflow. 

The method used to determine whether an x^ is in the 
regression at this step Involves ascertaining the sign of V^. This 
is due to the mathematical method of elimination which is used on 
the symmetric matrix. 

If is negative, 3^ is in the regression at this stage. 
The smallest value of is found for all in the regression. 
Denote this as and the corresponding subscript i as Nmi«. If 
is positive, is not in the regression yet. The largest value of 
is found. Denote this as '^max corresponding subscript i 

^max* 
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The following criteria are used in selecting the variable to 
add to or remove from the regression. 

(1) If ( |vj^ I • ^ the corresponding variable Xj^ 

is insignificant at this time, therefore will be removed 
from the regression before adding any other variable and k-|^ 




min** 


(2) If OTniax * 4^/ (» "V .)>F1, the corresponding variable 

'■n,n max^ 

x^ is significant and will be added to the regression. We 


then choose k 
The pivot a, 


N_ 


will be chosen according to the above rules. 


kl,kl 

FI and F2 are F values which should be preselected for a 
problem. These values are based on the approximate con- 
fidence level. The program will run into a loop if these 
values are not chosen correctly. 

2«3*5*b* Matrix Operations 

Given an element a which is the pivot, selected according 
k,k 

to the criteria in 2.3.5*a> where k corresponds to the 
subscript of one of the independent variables. The new 
matrix will be generated by the following algorithm. Denote 


the new matrix elements as a 




The new 





for i-j-k 

(12) 

*kj/*kk 

for i-k, j/k 

(13) 

®ij ” ^k *kj /*kk 

for i/k, j/k 

(lU) 

“®ik/\k 

for i/k, j-k 

(15) 
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2,3.6 Calculation of the coefficients of regression equation 


After each variable entered or removed from the regression, the 
following are computed: 

The standard error of the dependent variable: 


a / i> 

nn 


\] 

The regression coefficients: 


‘»i • «in • n/cTi 
The constant term; 

^0 “ y ’I^i *i 

The standard deviation of the regression coefficients: 


«b. 


\ ^ii/n 


where i corresponds to the subscript of the variable in the 
regression at this stage and n the subscript of the dependent 
variable, is the element of the matrix at each stage. 


( 16 ) 

(17) 

( 18 ) 

(19) 


2,3.7 Calculation of predicted values and deviations 

The computation of the predicted values for dependent variables is 

based on the final regression equation. The deviation between the 

actual and predicted values of the dependent variable is 

for k - 1, 2, ... m 
k 


3, Minimum System Components Reqvdred 
3*1 Fortran Compiler D 
3*2 One work tape. 


4. Input 

The input data Is read from the card reader. The input data 
format will be given in detail in a later section. 

5 . Output 

All output is printed on-line. 

6. Subroutine 

The square root routine must be supplied from the library of 
subroutines. 

7 . High Speed Memory Requirements 

Memory Space Notation Used In Remarks 

Program 

2 

A N VECTOR (I,J) Matrix A space has multiple 

use as; 

1) sums of squares and 
cross products 

2) residual sums of squares 
and cross products 

3) simple correlation 
coefficients 

4) partial correlation 
coefficients 

5) an inverse matrix 


1 ^ik 

N 

SUMX(I) 

(fl 

N 

SIGMA(I) 

^Ik 

N 

X(I) 

5i 

N 

XBAR(I) 

^1 

Sbl 

N 

N 

BETA(I) 

SB(I) 

DEV 

N 

STDEV(I) 


4-8 


7.1 Estimate of data size 

(N^+?N)(d+2) ^ M - 16,500 
where M Is the memory size 


d Is the digits of accuracy 


As an example; Assame d = 10 


Memory Size 
20K 
24K 
23K 
32K. 


Data Size 
< 14 
N < 22 
N ^ 28 
N < 33 


* Remark; This version is written for 16K machine with 
precision d = 10. In order to adapt this to 
other size machines, the arrays in the "common” 
statements of each chain must be changed 
according to the above estimated formula. 


8. Preparation of Input Data 


8.1 Size of numbers 

The largest sum of squares must be less than 10^. The size 
of numbers depends on the number of data points. The 
number that Is punched on a data card should not be more 
than 10 characters, including sign and decimal point. 

8.2 Number of variables 

This program will accept a maxlmuim of N Independent variables 
calculated according to the estimated formula in section 7 
and one dependent variable of a homogeneous set of data. 

8.3 Number of problems 

A maximum of 99998 problems can be run at one time. 

8.4 Number of data points 

The number of data points (observations) that can be 
handled is less than 10° and is limited by the restriction 
that the largest sum of squares must be less than 109, This 
program handled an equal number of homogeneous data points 
for each variable. 
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8.5 Decimal point 

The decimal point of input data need not be entered if 
the immber is an integer. The decimal point i-iill be 
assumed to be at the extreme right. For example, columns 
21 to 30 are used for the variable Xt_, If a ”25'* is 
punched in columns 28 and 29, then xj will be 250. A 
number vjith a decimal point can be punched in anywhere 
viithin the specified columns, for that variable. For 
example, a 25. in column 21 to 23 is the same as if 
punched in 28 to 30. 

8.6 Card format 


3.6,1 Control card of each problem 


Columns 

Contents 

Description or Remark 

1-5 


Not used 

6-10 

Problem ID, 

An integer < 99999 

11-15 

Number of Obser- 
vations 

An integer 

16-20 

Number of vari- 
ables 

An integer (see 7 
about its size limit- 
ation) 

21-30 

FI 

Decimal 

31-40 

F2 

Decimal 

41-50 

TOL 

Decimal 

51-BO 


Not used 

Data card 



Columns 

Contents 

Description or Remark 

1-5 

Blank 


6-10 

Input Data ID 

This no, should be 


exactly the same as 
the no, which appears 
in Column 6-10 in the 
Control Card, 
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Columns Contants Description or Remark 

11-15 Logical Sequence An integer. 1 stands for 

Number first 5 variables data 

card, 2 for second five 
variables data card, etc; 
if number of variables is 
less than 5» then a 1 is 
sufficient for all data 
cards. 


16-20 Observation No. An integer. 1 stands for 

first observation, 2 for 
second observation, etc. 


21-30 

^1 

Decimal, 

value 

of 

first data 

31-'^I0 

^2 

Decimal, 

value 

of 

second data 

41-50 

X3 

Decimal, 

value 

of 

third data 

51-60 

X4 

Decimal, 

value 

of 

fourth data 

61-70 

^^5 

Decimal, 

value 

of 

fifth data 

71-80 


Not used 





8.6.3 Last control card 

This card is needed at the end of a set of problems. 


Columns Contents 
1-5 

6-10 The Integer 
99999 

11-30 


Description or Remark 
Not used 


Not used 
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8.6. U Card Sequence 


2nd Probleirt 
data 


First Problem 
data 


last control 
card 




9. Sample Problem for Step-wise Regression 


9.1 Problem 2 

Find the best fit of 

y - . bi ♦ bj X2 ♦ bj Xj ♦ b|^ x^ 

with the following avai-lable data of 13 observations of values of 


5 variables: 




*2 

X3 

XU 

Y 

1 

7. 

26. 

6. 

60. 

78.5 

2 

1. 

29. 

15 . 

52 . 

7U.5 

3 

11. 

56. 

8. 

20. 

10 U .3 

li 

!!• 

31. 

8. 

U7. 

87.6 

5 

7. 

52. 

6. 

33. 

95.9 

6 

11. 

55. 

9. 

22. 

109.2 

7 

3. 

71, 

17. 

6. 

102.7 

8 

1. 

31. 

22. 

lUi. 

72.5 

9 

2. 

5U. 

18. 

22. 

93.1 

10 

21. 

U7. 

U. 

26. 

115.9 

11 

1. 

UO. 

23 . 

3U. 

83.8 

12 

11. 

66. 

9. 

12. 

113.3 

13 

10. 

68. 

8. 

12. 

109. U 


9.2 Input Data Cards 



fiOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOIIQOOOOOOQOOOOOOOOOOOOOOOOODODOOODO 

M 3 4 S ( 7 I S 10 II n » 14 19 l( i; II n»?l 22 2324 7S»n 212130 3133 33 34 3S34 J7M39 40 41 42 43 44 4S«< 47 4|49»l$>M)i»«»$(579l99(«tinUMBK(7»t»30n 7233 74 71 7{n 71 7110 

M 1 1 11 n 11 1 t|i 1 1 1 1 11 1 1 1 n n 1 1 1 M 1 11 11 n n n I II i|n f 1 1 1 1 n 1 1 1 1 1 1 ti 1 1 11 1 1 1 1 1 1 1 1 n 1 

22222222212222222 222222212222222222222222222222222222222222222222227222222222222 
3 3 3 3 3 3 3 3 3 3 3 3 3 1 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 1 3 3 3 3 3 3 3 3 3 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 
44444444444444444444444444444444444444444444444444444444444444444444444444444444 
5S555555555555S5555|55555|5S55SS555|5S555S5S5S555555555S5555!i55S555S55555555555S 
66B686866686666S66S6G6S66666666666664666G66fi6S66E66G66fi666666e6G66B8£8BE6G666E6B 
77777777777777777777777777777777777777777777777777777777777777777777777777777777 


999999999999S99999999999999999999S999999999999g9393999999S999d9S9S39599999599999 

1 2 3 4 S 1 ; I 9 ID 11 12 13 14 1$ 11 17 19 20 21 22 2} 24 29 21 27 21 21 30 31 32 U14 2SX }7 31 39 40 41 42 *3 44 494( 47 41 49 90 i< 92 93 94 99 U 97 91 91 H (I «2 (3 H 69 (((?('> *9 70 71 12 73 74 K 70 77 79 79 10 
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9,3 Output Data 


STEP-WISE REGRESSION ANALYSIS PROGR AM ( 5TR AP I ) 


PROBLEM NO, 2 

NO, oP OBSERVATIONS = 13 

NO, OF variables s 5 

F LEVEL TO ENTER VARIABLE = 2,500 

F LEVEL TO REMOVE VARIABLE » 2,500 


SUM OF variables 

X( 1) s 97,000 XI 2) = 626,000 X( 3) » 153,000 

XI A) I 390,000 

Y s 1240.499 


RAW SUMS OF SQUARES AND CROSS PRODUCTS 


X ( 

n 

VS 

X I 

1) 

8 

1139,000000 

X I 

i> 

VS 

X( 

2) 

s 

4922,040000 

X ( 

1) 

VS 

X 1 3) 

S 

759,000000 

X ( 

1) 

VS 

X ( 

4) 

8 

2620,000000 

X ( 

2> 

VS 

X ( 

2) 

s 

33050,000000 

X 1 

2) 

vs 

X 1 3) 

8 

7201,000000 

X ( 

2) 

VS 

X I 

4) 

8 

15739,000000 

X I 

3» 

VS 

X ( 

3J 

= 

2293,000000 

XI 

3) 

VS 

XI 4) 

S 

4628,000000 

X ( 

4) 

VS 

X ( 

4) 

8 

15062.000000 

X 1 













XI 

1) 

vs 


Y 

8 

10031.999940 

X| 

2) 

VS 


Y 

s 

62027,799970 

XI 

3) 

vs 

Y 

s 

13981,499940 

X I 

A) 

vs 


Y 

8 

34733.299970 















Y 

vs 


Y 

8 

121086. OB9900 
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X( 1) 
U 

r 


7*461 

30,000 

93,423 


average value or variables 

X( 2) • 40*133 


X( 3) 


11*769 







RESIDUAL 

SUMS 

or 

SQUARES 

AND CROSS PRODUCTS 








X ( 

n 

vs 

XI 1) 

• 

413.230766 

X( 

u 

vs 

X( 2) 

m 

291.076922 

XI 

11 

VS 

XI 

3l 

8 

•372,619364 

XI 

1) 

vs 

XI 4) 

■ 

•289,999999 

X< 

21 

vs 

XI 21 

m 

2909,692305 

XI 

2) 

VS 

X 1 

3) 

B 

-166,938461 

X ( 

2) 

vs 

XI 4) 

B 

•3040,999993 

XI 

31 

vs 

XI 31 

m 

492,307692 

XI 

3} 

VS 

X 1 

4) 

S 

36.000000 

X I 

4) 

vs 

XI 4) 

B 

3362.000000 

XI 













XI 

1) 

vs 

Y 

B 

773,961336 

XI 

2) 

vs 

Y 

B 

2292.993644 

XI 

3} 

VS 


Y 

B 

-616.230760 

X 1 

4) 

vs 

Y 

B 

•2461,699996 














Y 

vs 

Y 

B 

2719,763076 















X( 1) 
X( 4) 
Y 


standard deviations or VARIABLES 

9,682 XI 2) > 13,360 X( 3) s 6.403 

16.736 
13*043 


CORRELATION COEmCIENTS 


X 1 

11 

VS 

XI 11 

B 

1.000000 

X 1 

11 

VS 

XI 41 

8 

•.249449 

XI 

21 

vs 

XI 41 

8 

•*972994 

X 1 

4) 

vs 

XI 41 

B 

1,000000 

XI 

1) 

vs 

Y 

8 

,730717 

XI 

4) 

vs 

V 

B 

•.621309 


Y 

vs 

Y 

a 

1.000000 


XI 

n 

VS 

XI 21 

a 

,226579 

* 1 

21 

VS 

XI 2) 

B 

1,000000 

* 1 

31 

VS 

XI 31 

B 

1,000000 

X 1 






X 1 

21 

VS 

Y 

B 

,616252 


XI 

1) 

VS 

XI 3) 

B 

•.824133 

XI 

21 

vs 

XI 3) 

9 

-.139242 

XI 

3) 

vs 

XI 4) 

9 

.029537 

XI 

3) 

vs 

Y 

a 

-.934670 


standard error Or V « 13,043722 



STCP HO, 

1 







variable 

ENTERING 

X- 4 






F level « 

24, 

87U 






STANDARD 

error OF 

Y « 6,9639 






corstaht 

TENH ■ 

117,56793 






variable 

NO, 

COEFFICIENT 

STD 

ERR 

OF 

COCFF 


X- A 


•,73616 




,15499 


step no. 

2 







variable 

entering 

Xm 1 






F level 1 

119, 

0462 






standard 

error of 

Y ■ 2,7342 





H* 

>4 

CONSTANT 

TERM • 

103,09738 






variable 

NO, 

coefficient 

STO 

ERR 

OF 

coeff 


X- 1 


1,43995 




.13841 


X- 4 


•.61395 




,04664 


STEP NO, 

3 







variable 

entering 

X. 2 






F level « 

9, 

5842 






standard 

ERROR OF 

Y > 2,3087 






CONSTANT 

TERM • 

71,64830 






variable 

NO, 

COEFFICIENT 

STD 

ERR 

OF 

COEFF 


X- 1 


1,49193 




,11699 


X- 2 


,41610 




, 18561 


X- 4 


•,23694 




,17328 


STEP NO. 4 

VARIABLE REMOVED X* 4 

r LEVEL • U6632 

STANDARD ERROR OF Y ■ 2.40S3 

constant TERN • 52.5TT34 

variable NOi coefficient STD ERR OF COEFF 

X- 1 1,46830 ,12130 

X- 2 ,6622S ,04S6S 


predicted v5 actual results 


RUN 


NO* 

ACTUAL 

predicted 

deviation 

1 

T8. 50000 

80,OT400 

-1*57400 

2 

T4, 29999 

T3. 25091 

1*04908 

3 

104,29999 

105.814T3 

-1*51473 

4 

8T, 59999 

89*25847 

-1 *65847 

5 

95,89999 

97,29251 

-1*39251 

6 

109,19999 

105.15248 

4,04751 

T 

102.69999 

104,00205 

-1.30205 

8 

T2, 50000 

74,57541 

-2,07541 

9 

93,09999 

91.27548 

1,82451 

10 

115,89999 

114,53754 

1.36245 

11 

83,T9999 

80.53567 

3,26432 

12 

113.29999 

112,43724 

.86275 

13 

109,39999 

112,29343 

-2,89343 


END OF calculation* PPOBLEM NO* 2 

END olympiad 


9.4 Special Comments on Error Conditions 

If the following comments appear on output: 

9.5.1 data id does not agree with control card. 

INID 4 XXXXX XXXX NOC = XXXXX INKO = XXXXX 

9.5.2 DATA CARD OUT OF LOGICAL SE:^UENCE. CHECK. 
INID = XXXXX XXXXX NOC = XXXXX INNO = XXXXX 

9.5.3 DATA CARD OUT OF NUMERICAL SEQUENCE. INID ^ 
XXXX XXXXX NOC = XXXXX INNO = XXXXX 


then these errors are caused hy the mispunching of the 
identification number of the test data of a problem or 
misplacing of data cards. The program will continue, 
but the user should recheck the input data deck. The 
nvimber INID, NOC and INNO will help the programmer to 
pinpoint the mispunched card or misplaced card. After 
correcting the input data, the program should be re-run. 
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SECTION V 


RANDOM NUMBER GENERATOR D PROGRAM 


I. Introduction 

1.1 Purpose 

To generate a set of random numbers, 

1.2 Method 

1.2.1 Symbol 

n = nvimber of random numbers to be calculated. 

= jth random number. 

1.2.2 Description 

The routine uses an additive process: 

The aj , J=l, ..., 16 are specified as input data. 
These 16 random numbers may be obtained from any 
table of random numbers (e.g. Quality Control and 
Industrial Statistics , by A. J. Duncan). After a 
set of random numbers has been generated by this 
routine, the last 16 numbers obtained may be used as 
input to a future run, to generate a new set of 
random numbers. 

II. Usage as a Routine 

2.1 Minimum System Components Required 

2.1.1 16k fortran 

2.1.2 One on-line card reader 

2.1.3 One on-line printer (if printed output is desired) 
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2.1.4 One on-line card punch (if punched output Is desired) 

2.1.5 One tapet Program Tape 


2.2 Input 

The Input data Is read from the card reader. The first card 
must be a control card, followed by four data cards. 

2.2.1 Control Card Format 
Columns Contents 


Blank or «1" 


Description or Remark 

Not used. 

Blank for printed 
output, "1" for 
punched output. 

Not used. 

An integer, right 
Justified, specifying 
the number of random 
ntimbers to be gener- 
ated. 

56-80 Not used. 

2.2.2 Data Card Format 

Columns 1-25 26-35 36-40 41-50 51-55 56-65 66-70 71-80 


1-30 

40 

41-45 

46-55 


N 


1st data card 

ai 

®2 

^3 

*4 

2nd data card 

»5 

®6 

&? 

*8 

3rd data card 


aio 

*11 

*12 

4th data card 

^13 

»14 

*15 

*16 

The aj , J«1 , . . . , 

16 are 

10 decimal digit 

random 

numbers. ' 


decimal point is assumed at the extreme left. The other 
colimns are not used. 


2.3 Output 

2 . 3.1 Printed 

The generated random numbers will be printed on-line 
If specified on the control card. The output will 
Include the first 16 random numbers, used as Input. 

2 . 3.2 Punched 

If punched output Is specified, then the new random 
numbers (excluding the I 6 Input random numbers) will 
be punched on-line onto cards. The format of these 
cards will be Identical to the format of the Input 
data cards (See 2.2.2). 

2.4 Sample Input Data 



oooooooooooo(ioooooDooooooooooooooooooooooooooouoi]Qoooooooooo''oooooaoaooo(iooooooD 

• ? 3 4 S 6 7 B 9 10 n 13 13 IB IS 16 17 IB IB 70 31 27 73 24 2S2E 77 28 29 30 37 33 33 34 353E 3 ; 3S 39 40 41 42 43 44 4S4E 47 48 49 S0S1 53 53 S4S5SG S7 58 S96QG1 E7fi3 64 ESE6 67 6B 69 7f7 7l 72 73 7B7S76 n 7B79BB 

M 1 1 1 n n 1 1 1 1 n 1 11 111 11 1 1 11 11 11 1 ' 1 1 1 n 1 1 : 11 1 1 11 1 11 1 1 1 n 1 1 n 1 n 11 1 1 1 1 n n 1 1 1 1 n 1 1 

2222222222222222222222222222222222222222222222222222222222222222222J22222222222? 
33333333333333333333333333333333333333333333333333333333333333333333323333333323 
M44444444444444444444444444444444444444444444444444|444444444444444444444444443 
3S55S5S5S5555555955555555S5553555SS555555555S5S555SS55S55SS3555555SS5555S553555b 
66G68G6E6E6E6666E6E666E6EB6E66E666666E866SE66666E66666666e6666668E6e6666666EEB66 
7717777777777777777777777777777 7777777777 7777777?? 77777 7 7 777 77777 777 7777777777 77 
8 8 8 3 8 8 8 8 8 8 8 8 i 8 8 8 3 S 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 S 8 8 8 8 8 8 8 8 8 1 8 8 8 8 8 8 8 3 8 8 8 8 8 8 8 8 8 8 8 e 8 8 8 8 8 
33939999999999999959S9999989929S999SS99999‘i99S993999S|939933ii999S3999999933S999S 

• 7 3 4 5 G ’ 6 3 ID III? 13 74 15 15 me l» m <1 ?2 23 24 25 2C .'7 28 29 3C 3» 3? 33 34 3S 36 37 38 3? 40 41 «2 <.3 44 45 45 47 dj 50 5i 5’ «7 54 55 56 57 5S 5S 50 5' 6? 63 f : 63 56 M 58 69 70 71 72 73 74 75 T? T7 78 TJ 60 
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The 16 random numbers above were taken from Quality Control 
and Industrial Statistics , by A. J. Dunoan. The card numbers 


in column 20 are inserted merely for convenience. 

Usage as a Subroutine 

3.1 Minimum System Requirements 

3.1.1 16k fortran 

3.1.2 One tapet Program Tape 

3 . 2 Call Statement 

CALL RANDOM (J, N, L, M) , where J, N, L, M are a list of 
actual arguments. 

J must be an integer constant, J«3 

N is an integer constant, specifying the total number 
of random nmbers to be generated. 

L must be an array of integer variables of dimension 8, 
containing random numbers (ai, a.2* ^8^* 

M must be an array of integer variables of dimension 8, 
containing 8 distinct random numbers (a^, a^Q, ..., 

*16 )• 

3 . 3 Output 

With each call to the subroutine, RANDOM will deliver 16 
random numbers. These will replace the previous 16, stored 
in L and M. N will be decremented by 16, after each call. 
Thus, when N 5 the specified number of random numbers 
will have been generated. 


IV, Disclaimer 

This program has been checked but no guarantee Is expressed or 

Implied to the user concerning Its functioning, Honeywell will 

✓ 

cooperate In making functional those programs written by them. 
However, they assume no responsibility for programs distributed 
but not written by Honeywell, 
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